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Abstract 


The movement of ionic solutions is an essential part of biology and technology. Fluidics, from nano- to 
micro- to microfluidics, is a burgeoning area of technology which is all about the movement of ionic 
solutions, on various scales. Many cells, tissues, and organs of animals and plants depend on osmosis, as 
the movement of fluids is called in biology. Indeed, the movement of fluids through channel proteins (that 
have a hole down their middle) is fluidics on an atomic scale. Ionic fluids are complex fluids, with energy 
stored in many ways. Ionic fluids flow driven by gradients of concentration, chemical and electrical 
potential, and hydrostatic pressure. Each flow is classically described by its own field theory, independent 
of the others, but of course, in reality every gradient drives every kind of flow to a varying extent. 
Combining field equations is tricky and so the theory of complex fluids derives the equations, rather than 
assumes their interactions. When field equations are derived, rather than assumed, their variables are 
consistent. That is to say all variables satisfy all equations under all conditions with one set of parameters. 
Here we treat a classical osmotic cell in this spirit, using a sharp interface method to derive interface 
boundary conditions consistent with all flows and fields. We allow volume to change with concentration, 
since changes of volume are a property of ionic solutions known to all who make them in the laboratory. 
We consider flexible and inflexible membranes. We show how to combine the energetics of the membrane 
with the energetics of the surrounding complex fluids. The results seem general but need application to 
specific situations of technological, biological and experimental importance before the consequences of 
consistency can be understood. 


1 



Introduction 


Osmosis moves ionic solutions throughout biology and technology. It is hard to find a more 
widespread phenomenon. All biology occurs in ionic solutions 1,2 that move on many scales, including 
atomic scales smaller than nanometers, and a great deal of chemistry involves ionic movement as well. 
The modem technology of fluidics—macro, micro and nano—usually moves water and ions and thus 
involves osmosis at every scale. Indeed, modern names “ * fluidics ” may be viewed as a renaming of 
a classical, 3 if not ancient word 4 into more modern language. 

Ionic solutions involve energy stored in many forms, pressure, concentration, electric and 
electrochemical potential, steric interactions, and chemical energy. 1 Ionic solutions are idealized as 
perfect gases 5 in elementary textbooks, whether biological, 6 technological (polarography) 7 - 8 or 
general. 9 ' 16 Ionic solutions idealized this way do not have stored energy. The idealizations also do not 
involve time (in most cases) so they cannot deal with flow or friction. 

Ionic solutions are complex fluids, 17 " 19 not simple at all. 

Classical analysis of ionic solutions started with a thermodynamic and thus equilibrium analysis 
without flows. 11 ' 16 Thermodynamics was dynamic in its historical importance as a turning point in the 
history of science 20 but it was not dynamic—it did not involve time— in its physics or mathematics. 
Without time, it could not include flow or friction, at all. Flows occur in almost all living systems, and 
technological devices, so thermodynamics was extended ad hoc to deal with flows, 21 " 24 often with great 

25-30 

success. 

The theory of complex fluids is a dynamic theory that always includes flows of many types, 
usually in a coherent consistent way. (‘Consistent’ means that all variables satisfy all equations and all 
boundary conditions with one set of parameters). In complex fluids like electrolytes, ions interact and 
move as components of complex fluids. 18 ' 3137 The solution itself flows more or less as water itself 
would move (without the ions). Water molecules and ions move (partly) by bulk flow, that is to say, 
they move (partly) by convection described classically by Navier Stokes equations. Water molecules 
and ions also move (partly) by diffusion. Ions move (partly) because of the electric field. Water also 
moves in an electric field because of dielectrophoresis 38 ’ 39 when the electric field is nonuniform i.e., 
\VE\ > 0, where E is the electric field. Each of these flows varies with location and is described by 
field equations—typically partial differential equations in space and time—along with boundary 
conditions, that idealize the physics of the particular setup in which the flows occurs. One example of 
these field equations are the Maxwell equations of electrodynamics. 40 ' 42 Few are more important 41 - 43 . 
Most of our technology 44 ' 49 is built from the Maxwell equations, perhaps because they are so precise 
and universal. 

But even Maxwell’s equations need to be extended to deal with osmosis. The Maxwell’s equations 
do not involve hydrostatic pressure or concentrations of chemical species, for example, and so they give 
the same residts, no matter what the pressure or concentration gradient. Thus, Maxwell’s equations in 
themselves cannot deal with many crucial phenomena in biology and technology. 


1 We define ‘chemical’ forces between atoms as those that significantly change the spatial distribution of electrons in the 
atoms. In this definition, dielectric interactions are not a chemical force. More precisely, force produced by the induced 
charge, that is proportional to the local electric field, is not a chemical force. Dispersion forces arise from the quantum 
fluctuations in induced (‘dielectric’) forces and so are not considered chemical in this definition. 
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Pressure and concentration are responsible for many phenomena, in the laboratory, in technology 
and in biology. 21 24 Experiments show that gradients of hydrostatic pressure and concentration change 
flows a great deal, quantitatively and qualitatively in physical and living systems. 11 - 21-24 - 50 Maxwell’s 
field equations cannot change with changes in hydrostatic pressure or flow because the field equations 
do not contain variables describing pressure or flow. Classical electrodynamics must be extended if it 
is to deal with the range of forces and flows so important in technology and life. 

Conservation laws easily define the field equations of osmosis when a single force is involved, 
with some help from constitutive laws and boundary conditions (that describe the experimental or 
natural setup, more than anything else). 51 ' 52 But the situation is more complex when diffusion, and 
electrical migration, and bulk flow are all involved. 5 ’ Bulk flow moves ions and water. Diffusion 
changes electrical potentials, electrical migration changes diffusion, and potential and diffusion (even 
of nonelectrolytes and water itself) change and are changed by bulk flow. Technology and biology use 
ionic solutions in devices where bulk flow, diffusion, and electrical migration are important. A theory 
must then deal with most of these conditions, robustly with one set of parameters, if it is to be useful. 
Only robust theories can create technologies as powerful as our semiconductor electronics. We reach 
to that success 54 " 56 as our goal. 

Coupled field equations that satisfy all conservation laws with one set of parameters are consistent 
equations. Writing consistent equations is difficult. Even choosing the appropriate variables causes 
problems: for example, the operational definition of appropriate variables for an infinitely dilute 
sodium chloride solution remains requires extensive discussion 57 without definite resolution (after 
nearly a century of work). 

Field equations have different form when different variables are chosen. Equations are typically 
written in terms of the mass of each specific chemical species, for example in terms of the mass of 
sodium and the mass of chloride in a salt solution. But a linear transformation would allow the 
equations to be written in terms of total mass density and net charge density, instead of sodium and 
chloride concentrations. The equations, coupling terms, and boundary conditions would be different. 
Choosing the ‘right’ set of equations is difficult. 

Experimental results do not depend on how we define variables or formulate theories. The goal of 
theoretical analysis should be robust transferrable theories—that fit a range of data over a range of 
conditions with one set of parameters. The semiconductor technology which has remade our world is 
catalyzed by its successful theories. Similar robust theories of osmosis are likely to catalyze 
applications in biology and with substantial effects. 

The theory of complex fluids has developed consistent theories for many complex systems 35 - 58-60 
that are mixtures (or even ionic solutions) with components that store energy. Variational methods 
deal successfully with magnetohydrodynamics systems 61 , liquid crystals, polymeric fluids 36, 37 , 
colloids and suspensions 58 ' 62 and electrorheological fluids 63 . Variational methods describe solid balls 
in liquids; deformable electrolyte droplets that fission and fuse 58 - 64 ; and suspensions of ellipsoids, 
including interfacial properties of these complex mixtures, such as surface tension and Marangoni 
effects of ‘oil on water’ and ‘tears of wine’. 58 ' 62 ’ 65 ' 66 

The Energetic Variational Approach (EnVarA) of Chun Liu, et. al., has been applied to ionic 
solutions of interest in references. 19 ’ 67-72 EnVarA combines the Least Action Principle of Hamiltonian 
dynamics with the Maximum Dissipation Principle of Onsager (who used the dissipation function of 
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Rayleigh) into a set of field equations, typically in Eulerian coordinates. Variations are taken with 
respect to two different variables, position and flow rate. The resulting field equations are written in 
one set of coordinates—usually Eulerian—with push-forward and pull-back methods. 

EnVarA is an approach, not a general method. An approach adapts to the experimental setup and 
questions of interest. It avoids the abstraction of a universal formulation that may prove hard to 
understand 73 ' 76 and hard to apply uniquely in specific cases. 

The theory of complex fluids requires models that describe the constraints put on the field 
equations by biology, technology or experiments, as well as the field equations themselves. It requires 
boundary conditions as well as field equations. The constraints and boundary conditions often are the 
dominant determinant of the dynamics of the whole system. Describing these constraints with 
mathematical boundary conditions can be challenging, as in the moving contact-line problem 77 79 and 
when describing classical setups of the theory of stochastic processes. 80 " 84 

Boundary (or interface boundary) constraints of the physical problem can be converted into 
mathematical boundary conditions using the energy variational method. 78 - 79 ' 85 " 87 The maximum 
dissipation principle (attributed to Onsager 88,89 20 " 22 developed from the work of Rayleigh 74 - 75 - 90 ’ 91 ) i s 
often used to describe sharp interfaces. 1819 Here we combine the energy variational approach (to field 
equations) with the sharp interface approach to physical constraints and use the combination to analyze 
a classical osmotic cell 3 - 51,92 —two baths separated by a semipermeable membrane—with flows driven 
by electrical, diffusion, and pressure fields, neglecting for the moment the steric and chemical forces 
of non-ideal solutions, or heat driven flows. Liu and colleagues have had some success creating what 
we call ‘thermal dynamics’ that deals consistently with heat driven flows in the spirit of EnVarA and 
the theory of complex fluids. 93 ' 95 

The variational approach is widely used to analyze complex fluids because it derives field equa¬ 
tions and boundary (or interface boundary) conditions, rather than assumes them. Once energy and 
dissipation functionals are defined by a physical model, the interaction terms in the field equations for 
each flow (bulk flow, diffusional flow, and electrical migration) are determined by algebra and analysis 
with minimal adjustable parameters. All variables satisfy all equations with one set of parameters when 
the analysis is done correctly. Results fit experiments, over a range of conditions with one set of pa¬ 
rameters, if the original functional equations—from which the field equations are derived—are a cor¬ 
rect enough model of the physics of the system, and the description of the constraints and boundary 
conditions is adequate. Results of the theory are then transferrable (using the language of chemis¬ 
try 96 - 97 ). Technology shows the importance of transferrable results. Our electronic/semiconductor tech¬ 
nology would hardly exist if devices had different properties in each location in the remarkably diverse 
circuits they use. 44 

It is important to understand that field equations and boundary (or interface) conditions written 
without variational methods are likely to be incomplete, truncated, and thus not transferrable. 
Mathematical models written without variational methods are likely to omit terms needed to satisfy all 
field equations and boundary (or interface) constraints with one set of parameters, particularly if they 
are written as equations with right hand sides equal to zero, as is often the custom. (Vanishing terms 
are easy to overlook when applying conservation principles!) When the omitted terms have important 
effects, the truncated field equations cannot satisfy all field equations with one set of parameters. They 


4 



are inconsistent. As the system is transferred from one set of conditions to another parameters must be 
changed to fit the data, and the new parameters are often only known by the data that is fit by them. 

Our analysis of osmotic flow is in the spirit of a variational approach, combining elements of 
EnVarA and a sharp interface analysis that itself uses variational ideas. We also deal explicitly with 
the variation of the density of a solution with concentration, 98 responsible for the difference between 
molal and molar concentration units, taught to everyone who makes up solutions (from solid solutes) 
in a laboratory. 11 - 99 This density dependence may be important in some experimental setups and 
conditions. 

Our approach needs to be applied to each experimental setup before one knows its importance: 
different setups will emphasize different terms. Many theories are found to apply only in special setups 
and a particular range of conditions. 7,8 100 102 Disagreements are found in the literature of osmotic 
flow—some are mentioned in the Historical Note of this paper. Perhaps some of these disagreements 
reflect the use of different inconsistent theories—which may omit different terms—to deal with 
different experimental setups and conditions. Many setups will be sensitive to density changes with 
concentration: molar vs. molal terms may be important quite generally. See our eq. (23)-(24). We are 
unaware of other consistent analyses of the density effect—that is so well known to chemists and 
biologists who have made solutions in their laboratories—although the literature is vast and so is our 
lack of knowledge of it all. 

A consistent theory should resolve differences, if the underlying physical model is correct. 2 
(For us the underlying physical model is written in functional equations and boundary (or interface) 
constraints, e.g., eq. (l)-(2) that follow.) We hasten to add that none of us are physical chemists, let 
alone experimental physical chemists, and we are certainly aware that there are many experimental 
issues involved we do not know much about and that need to be included in applications of our work. 

The paper is organized as follows. The mathematical model for a semi-permeable membrane is 
developed in Section 1. The development depends on several unknown variables. In Section 2, those 
variables are developed for the case where solution density is a function of ion concentration. Those 
unknown variables are evaluated for the case of constant density in Section 3. A historical section and 
discussion conclude the paper. 

Note to the reader . Variational approaches have not been as widely used as they might be, perhaps 
because they involve a large amount of tricky mathematics that has to be done and checked. We present 
all steps in our analysis, hoping to make the approach easier to use. Apologies are extended to those 
who find the detail overwhelming, or irritating. 


2 For us the underlying physical model is written in functional equations and boundary (or interface) constraints, 
e.g., eq. (l)-(2) below. 
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Mathematical Model 


We consider the dynamics of a fluid with ions passing through a seim-permeable membrane in 
the traditional setup of physical chemists and biologists 21 " 24,50,92,103 following 104,105 . We use the sharp 
interface model 77,78,85,86,106 to derive detailed specific conditions on the membrane. Let denote 
two compartments separated by a membrane T where Cl — Cl + U Cl~ . (see Fig. 1). D 4 ^ is the electric 
displacement vector field of classical electrodynamics and cf is the distribution of the of i th (i = 
1, ...,1V) species of ion in domain H 41 , respectively. p ± and u 41 are the density and velocity of 
solutions in left and right compartments, respectively. 



Figure 1; Schematic of a flexible membrane. Solid line 
is initial shape of membrane and the dashed line denotes 
the deformed membrane, n is the normal vector of 
membrane T from the left compartment Cl + to the 
right compartment Cl ~. 


Based on the Maxwell’s equations and conservation law for the mass of each ionic species, we have 
the following equations 


D ± = ^ ztecf, 


dct 


dt 


dp 4 

dt 


l - + V- (cfu±) = -V • j± 


in Cl¬ 


in fF 


( 1 ) 


'du 4 
. dt 


+ (u 41 • V^u- ) = V • (cr^ + crl), in H 4 


+ V • (p ± u ± ) = 0, 


in H 4 


Zj is the valence of the i th species, e is the elementary charge, j f is the flux of the i th species 
<r^ is the viscous stress, and <j| is the electric stress in domain Cl ± , respectively, and cr^ are 
symetric tensors. 

If we neglect magnetic forces, there exists an electric field E ± and electrical potential 0 ± 
such that the dielectric constant D 41 = — —£ 0 £ptj(f) ± i where £ 0 f r" is i n domain fl 41 , 

respectively. We assume that the dielectric constant £p (units: dimensionless) and the permittivity 
£ 0 £f- (units: farads/meter) are each a single real number. 
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We do not deal with the nonideal properties of ionic solutions, despite our understanding of 
their significance, and the importance of nonideality in general. 1 ’ 43 ' 1<)7 ’ 108 We have to start somewhere 
and we have our hands full with the mathematics needed to deal consistently with the idealized cases. 
When we do reach to include nonideality, we anticipate difficulties. Even the proper formulation of 
the field equations (with flow) in the nonideal case remains an open question.If we neglect the 
magnetic forces, there exist an electric field E ± and electric potential <p ± such that D ± = 
£ 0 spE ± — —EqE^V^, where £ 0 £p is the dielectric constant in domain fl*, respectively. We 
assume the dielectric constant £p (units: dimensionless) and the permittivity £ 0 £p (units: 
farads/meter) are each a single real number and do not deal with the nonideal properties of ionic 
solutions, despite our understanding of their significance, and the importance of nonideality in 
general 1 ’ 43 107 108 We have to start somewhere and we have our hands full with the mathematics needed 
to deal consistently with these idealized cases. When we do reach to include nonideality, we anticipate 
difficulties. Even the proper formulation of the field equations (with flow) in the nonideal case remains 
an open question. 

For the interface condition of the surface F, also based on conservation law 109 , we assume that 


c 



All on T < 


D± • n = Q{x), 

• n — p ± v ■ n — Q p , 

U~ • T — V ■ T, 

[(Tjj + <7 e ] • n — [Qptr] — F mb 


( 2 ) 


)f n + cf (u± - v) ■ n = Ji (x) 


where X(-,t) is the trajectories of points on the membrane, v — v T + v n n is the velocity of 
membrane with normal vector n and tangential vector x, [/] = f + — is the difference of / 
across the interface T. We also assume that = v — 0 on dT, i.e., we assume the edge (or 
boundary) of membrane is fixed. Here Q(x), Q p (x),and Ji(x) are surface charge density, solution 
flux and ion flux on the interface, respectively. F mb is the membrane force induced by mechanical 
and dielectric properties of membrane. 

In the next section, we will derive the explicit forms of j f, oj, Q (x). Q p , J L (x) and F mb 
based on the second law of thermodynamics generalized to deal with flows driven by different forces. 


Remark F mb is taken as a surface force for simplicity. However, it can include the other types of 
forces such as elastic force and/or dissipation force induced by friction inside of membrane. 
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Total Energy Functional . In our case, and many others, 58,78,85,86,106 the total energy functional is the 
sum of the kinetic energy part E kin (p,u), internal energy E int (ci, (p) and surface energy E r (<p,y 0 ). 
The internal energy is composed of the electrostatic energy part E es (<p ), and the Gibbs free energy of 
ions E ion (Ci). The total energy functional is the precise statement of our field model (that is to say 
our physical model without boundary (dQ) constraints or conditions). And the total energy functional 
is given by 


7 tot 


- E kin + E int + E r 

— Ekin E es + Ei on + Ef 


-XL 


( e kin + e es + e ion) dx + J ( e r + Yo)dS 


(3) 



\u ± \ 2 + -E ± -D ± + k B T 


X 


dx 


+ 




dS 


where c 0 is a characteristic ion density, C m is membrane capacitance, y 0 is the membrane surface 
tension. In the following, square brackets always denote the jumps across the interface. 

2+ fn± f ± dx = f n+ f + dx + f n _ f~dx is the sum of integration in left and right compartments 

weighted by the corresponding value of /. 


The dissipation functional is defined as 

| 2t7±|D±| 2 dx + ^ j A ± \V-u ± \ 2 dx + ^j^ j Dyef\Vp.f\ 2 dx 


+ 


j G 1 ([fi i ])dS + j G 2 (Q p )dS 


(4) 


Where D v — (Vu + (Vu) T )/2 is rate of strain, r/ ± and A ± are the two Lame constants 109 , p.y = 
jlf (cf , 0 ± ,P ± ) a nd Dy- are the chemical potential and diffusion coefficient of i th ion. The first three 
terms in dissipation functional are the dissipation induced by fluid friction, volume change and ion 
diffusion in the bulk region. The last two terms are the dissipation induced by irreversible osmosis .on 
the membrane. The forms of G x (x) > 0 and G 2 (x) > 0, for any x, will be discussed later. 



Derivation: Density Depends on Concentration 


We start with fluid density as a function of ion density p = p(c±, c ••• c^). Based on the biological 
applications, here we consider the case that p weakly depends on the the ion concentration. 


Permeable and Flexible Membranes 

In this section, we consider the membrane is deformable (see Fig. 1). We first present a 
generalized Reynolds transfer formula 110 - 111 when the membrane is permeable 


T t ZL ( / d *=ZL 


n±(t) 


df d 

dt 


dx + 


/ [/> 


ndS 


(5) 


= ZL (^ +f±vu± ) dx -1 @ Q > ds ~l/ ±u±n ds ■ 

where L(tP ±dS =LnS*f* dS + LnB-f~ dS - 

Based on the results in, 77 ■ 112 - 113 we have the following formula for surface defined function /, 

f fdS = f (/ + fV r ■ v)dS — [ (?L + vVf + fV r -v)d 

Jrtt Jrfu Jrrt 'i ' 


d 

dt. 


S iiit + Vndn f + Vt ' Vr f + f Vr ' v ) 


f % + v n d n fdS + f fv T ■ n m dl - [ fHn ■ vdS 
d r(t)Jdr(t) d r(t) 


= f (^-£) dS- f fHn ■ vdS 

Jr(t) ' dt J J r(t) 


( 6 ) 


where V r = V — n(n • V) is the surface gradient operator, / = ^ + v ■ V/ is the material derivative 

d n f d f 

on the membrane, = — + v n d n f is the normal time derivative, n m is the unit outward normal 

vector of membrane at the edge dr, i.e. n m ■ n — 0, n m ■ dl = 0 and H — —V • n is the mean 
curve of membrane T. Here we used the fact that v = 0 on dT in our model. 

The following useful lemma states that the normal time derivative of normal vector on the surface 
only depends on normal component of membrane velocity. 
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Lemma Let n and v = v n n + v T be the outward normal vector and velocity of membrane F, 
respectively. Then we have the following result 

d n n 

= -v r v„. 


dt 


The proof of Lemma is presented in Appendix A. 

Now we begin to derive the full model. If we take the derivative of the energy (3), 

dE tot dE kin dE es dE ion dE r 

~ km + - - + + -r^ = h + I 2 + I 3 + / 4 


dt dt dt dt dt 


(7) 


For the first term I ± , by using the last two equations in Eq. (1), third interface condtion in (2) 
and Eq. (5), we have 


h= iil.!j p± ' u± ' 2dx 

|u ± | 2 dp ± 


-V, 


±2 dt 


dx + 




f 

[M 2 1 

1 

2 


+ / I -p ± |u ± | 2 V • u ± dx — I e^u^-ndS 
^ Jn± 2 Jdn 

— — ^ j (Tjj-. Vu^dx — 'y' j (Tj\Vu ± dx + j [n ■ (cr^ + <r e ) • u]dS 


[ 

[M 2 1 

Jr 

2 


QpdS + 7 lfci 


— j <7^:Vu ± dx — ^ J cr^:Vu ± dx + J [(n ■ + a e ) — Q p u) ■ u]dS 


( 8 ) 


+ 1 y QpdS j^p±^^(V-jt + c ^V-u ± )-V-u ±) jdx + I lb , 


where we used the notation d i p ± — (c±), the decompostion u = u n n + u T , u T is continuous on 


,+ _ d P 

Liuiaiiuii uif. 

the membrane and 

kb = f an (0± + ot) ■ n) • u±dS - f gn e kin u ± ■ ndS. 
The colon : represents the double dot of two tensors A: B = Y,p ^ij B LJ . 
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A Lagrange multiplier p c is introduced to assure mass conservation, when solution density is a 
function of ion density. 

The second and third terms in Eq. (7) are , with the help of first two equations in (1) and last 
two equations in (2), given by 



= X j n± X ^ ( _V ' ** ~ V ' W M± )) dx 



dD ± 

dt 




E &ion]v ' W dS 


■I LI V//f -Jt dx + I LI Cf (V/rf) • u— dx 


J n± 

+ 11 i 


+ - n± ‘ i 


~f ^HiUi-n + ciU-ri) 

- j k ’ n ] ds + / ^ es + eion ^ v ' nds + Izib 

■ILl V/if ■ jfdx + ILl Cf (V/rp) • u- dx 


J n± ‘ . 
± sz i 


J n± 4 . 
± sz i 


- ds ~ f 


dD 

*w n 


dS + 


j (&es E &ion ^^ di^i 


(9) 


v ■ ndS + I 21 b 


where pf = z i e<p ± + k B T (in 


9 -) 


and 
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l 21b — 



-njdS — J k-i(jf + c f u± ) ‘ n dS. 


Due to the first equation in (1) and second equation in (2), the second term in Eq. (9) is rewritten 

r 


ZLZ C '- 7 ^ U± 


± dx = 


Jn±*- 

+ n i 


I LI ZiecfV^ ■ u ± dx + I LI k B T\7cf ■ u— dx 


+ “ i 


+ “ i 


-IL • u-^dx — yf y kgTcf 1 V • dx 


V ■'fl±V 


+ J y? k B TndS + J y? k B Tc L u ■ n 


dS 


- ^ | D± (Vui) • V0±dx - ^ | D± (VV^i) •u±dx + J [QVcp ■ u] dS 

+ 1 (D 1 ■n)(V(f) ± ■u ± )dS 
Jdn 

-I LI k B Tcf V • u-dx + LI ksTcf^ ■ ndS + II k B Tci u ■ n 


± - n± ‘ i 


dS 


= Y [ e o e±(V0± (8) V0±):Vu ± dx + Y [ ^^(VlVcp 1 ! 2 ) ■ u± dx 

V V 2 

+ | [QV0 • u]dS + | (D± •n)^ 1 • u± ) dS 
J r 

-11,2 k B Tcf V • u—dx + LI k B Tcf^ ■ ndS + II? k B Tc t u ■ n 

|V0±| 2 


+ i 


dS 


= ^ J + e o e±(v0 ± ®V0 ± 


-/ : Vi^dx 


I LI 

± i 


k B Tcf V • u 1 * 1 dx 


+1 [e„« ■ n]dS + J [<2V0 ■ «]dS + IE k B TciU ■ n 

r r r L j 


dS + l22b> 


12 



where 


I 22b — | e es u± ‘ n dS + | (D^ • ri)(y<fi- ■ u- )dS + I y k B TcfU ■ ndS. 
Jda Jdn Jdn^-r 1 


( 10 ) 


By using the Lemma on the interface F(t), we have 


d n Q d n D d n n dD 

—— = —r~ ■ n + —— • D = — ■ n + v n d n D ■ n - V r v n ■ D, 
dt dt dt at 


which yields 


dD d n Q d n Q 

~dt' H = ~dT + VrVn ' ° ~ Vndn ° ' H = ~dT ~ e ° £rVrVn ' Vr< ^ ~ VndnD ’ n 

Then the fourth term in (9) could be rewritten as 


I 


dD 

0 aT n 


dS 


f d n Q f 

= J [0] -fa~ ds - J [ £ o^r^r^n ■ Vr0 + <pv n d n D ■ n] dS 
= j [0] ~^ ds + J l^r^r • (0V r 0)K dS - j [<fid n D ■ n\v n dS. 


( 11 ) 


Combining (9), (10) and (11) yields 


h + h ~ 


j V/tf -jfdx+ ^ J + e o e±fv0 ± ® V0 d 

I *T2 ^ _l /2 \ 


±12 


|V0±| 


-/ : Vu ± dx 


lil k B Tcf V • u~dx + J |Ye es + ^ k B Tc\ u ■ nj dS + j [QV0 • u]dS 


J n± * . 

± sz i 


J dS ~ J 


( 12 ) 


— J v r -(v r (^0 2 )) v-ndS + j [cj)d n D ■ n]v ■ ndS 

I 


J \_^es &ion * TldS + I 21 ? 

where I 2b = I 21b + I 22b . 

For the last term in Eq. (7), by using the surface Reynolds formula (6), we obtain 

d n 


h 


= Cm [0] -77 ([0]) dS- ((e r + Yo)Hn) ■ vdS. 

Jr ClL J r 


(13) 
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Combining Eqs.(8), (12) and (13), we have 
dE tot 


dt 


— — ^ j a^-.V^dx 


t~e 0 e± 170±(8)i70± 


|70±| 2 


/ : Vu^dx 


n\° 


+ X/ + X V ^’ ; ^ dX_ X | + Z k BT'cfV-u ± 


J[}± “ 
+ sz i 


+ n± i 


dx 


dx 


+ Z t^(l¥ ( "- y?+Ct±I7 ' M±) “ 7 - M± ) 

+ J [(n- (a v + <r e ) - Q p u) -u\dS + j (^e es + ^k B Tc^u ■ n 
+ | [QV0 • u]dS - J r ^i]/i ^ - j[cp]^-dS 

— j v r • (v r 0 2 )^jj V nd77 + / ' n ] v ‘ n 

-X p 


dS 



+ 


£es "b Gjon 

d n 


] v • ndS + I 

\\u n \ 2 ] 

/J Jr 

2 


QpdS 


j Cm ^ dt dS ~ / + 7o ^ V ' n d5 + + 7 2^ 

By using the fact that 


11 , ( rf l 7 ? p «)" ‘ 




• jfdx + 


Pc Iv) ;Vn 


dS 


(14) 
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+ l dn \^Zlk)*- ndS 


d iP\ s 4 


v-f 


—XL p 


dip 


JidS 


i Z Vc vi e o dS+ Li p ‘Z d -v) J ‘ ndS - 


diP 1 ) ,4 


Eq. (14) can be written as 
dE tot 

~=~) , 

>n± 


dt 


■I : Vu ± dx 


— — y I ai: Vu^dx 
T 1 Jii± 

_ z l + (^ _eoe K 70±070± ~^ 

+ X /+X ' i* dx ~ J Z ds 

± n± i r i 

Z 1 ± Z ^ ~ ?c± ( x " Z ^ ci )| 17 ■ u± dx 

J ^(“4 + < r e) ' n ~ Qp u + |e cs + ^ k B Tc, j n + QV<p\ ■ u 
£pic i + $d n D • n—V r • (v r 0 2 ) j 


+ 


+ 


dS 


I 


^es "b ^ton 


v -ndS 


diP l^nl 


l f X'dtf 

+ “Pc / — 2 "C£ + 

Jr A - 1 P 


-> P c 

i 

d n 


QpdS 


+ / ^ dt ^ Cm ^ -Q) ds ~ f ( e r + Yo)Hv -ndS + I b 

= ~ ^ J + a V : Vu±dx ~ Z / + ^' 7 0 ± ® 


I : Vu ± dx 


+ Z / + Z v/ ^' j^ dx ~J Z ^ ids 


+ 


z i ± (- z ^ - p ‘ ( j - z dj v' c 0) 17 ■ u± 


dx 


(15) 
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+ 


+ 


J | (tTjj + <7 e ) • n—Q p u + ^e es + ^k B Tc^jn + QVcj) 

J &e) ' ^ — QpU ” 1 ” ” 1 ” 'Y^k B Tc^Jn + QV0 

| e es + e ion - ^ \iiCi + <pd n D • n—V r • |v r (^y 1 0 2 ) j 


dtP , l^nl 


Q„d5 


+ i _pc Z^ ci+ 

i 

+ j [0] ^ (Cm[0] ~ Q) dS — J (e r + y 0 )Hv ■ ndS + I b 

IL* 


dx 


^ J IffZ-eoeffv#*® 70±-!^-L/j 1 : 7m± dx 


7|2± 


+ X /+X ■ ^ dx - j X ds 

± /2± t r t 

+ Z 1 ± (- Z ^ - p? ( x - Z ^ ci )) 17 - uJ 

+ J (Fmb + [QV 0 ]) • V T dS 

+ l 


dx 


((^77 + <r e ) • n - Q p u + (e es + Si k B T q)n + 


+ ( + 

'r 


• ni? n dS 


+ 


I 


^es ^ton 


e e s + ^^B7’c i jw + QV0 j 

^ Pt c t + 0d n D • n—V r • ^V r (yy0 2 ) j 


dtp , |w„| : 


+ i “ pc Z^ c ‘ + 

d n 


Q P dS 


r d n r 

+ [0] -77 (C m [0] — Q) dS — (e r + Yo)Hv n dS + 4 

Jr at J r 


■ v T dS 


■ nu n dS 


v - ndS 


nQp dS 


v n dS 
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-IL 


ai: dx 


e 0 ep 70 ± ® V cp z 


It r 


+ X / + X v ^' y ^ dx ~/X^-^ d5 

+ /2± i r ■' 


| 70±| 2 


/ : Fit* dx 




(16) 


+ X J + y X kBTCl± “ 

+ J ( F mb + [C?V0]) • V T dS 
+ J if mb ■ n + i F n ] - (e r + Yo)H)v n dS 

+ l \f\ QpdS + J r [0] di (Cm[0] “ Q) d5 “ 7 - 

where jlf = pf ~vi^> h = hb + hb + f dn (p£ Si ’ n 

F n = QVcP ■ n+ dik B T Ci + 2e es + e ion - Ci ^)-V r • (v r (^0 2 )) + P(d n D • n) 
= Qd n (p + (^2e es - ^ CjZ;e0 j -V r • |v r + p(d n D • n) 

= Qd n 0 + ^e o frl' 7 r0l 2 + e 0 e r \d n (p\ 2 - ^ c f z f e0 j 

-(f o frlv r 0l 2 + e 0 e r (pA r (p) + 0(5 n D • n) 


= - CjZje0 + e 0 € r <pA r cp - pd n D ■ nj 

— — (0V • D + e 0 e r (pA r (p — <pd n D ■ n ) 

= — (4>d n D ■ n + 0V r • D + e 0 e r <pA r <p — (pd n D ■ ri) 

— —( 0V r • D r — D ■ n pH + e 0 e r pA r p ) 

= QcpH 


(17) 


and 


F p = -Q p u n + n ■ (er v + o e )-n + ^e es + ^k B Tc^j + Qd n p - ^p c ^-Ci + p ^ 


(18) 


17 



Comparing with Eq. (4), by taking quadratic forms of G\ = ^ [/q] 2 and C 2 = we have 


j± = _BLfL V p± 

Jl k b t 


= 2t7±D± + A 1 (17 • i^)/ - pi 


'v 

r± — 


^ = SoSpV^ ® 170± -^|70±| 2 / 

9i 


• + g 0 g r 




5iP± 


P ± = Si k B Tcf +p± (l -Ef-^c, 

<3 = C m [0] 

Q p = -p 2 /f(*){[^]} 

F-mb ■ n — ( e r + 7o)tf - lA] 

r-F mb = -Q[V r 0] r 


f) 


in 

in 

in flA; 
in fi 1 

in 

on r 
on r 

on r 
on r 


(19) 


where /f(x) and g L are the permeability and conductance of the membrane, respectively. 

Remark Here we took G i and G 2 as the simple qudratic functions, i.e. G x = Y,i(di/e 2 )[P-i\ 2 and 
G 2 = Qf,/K(x). In fact, any flux J (x) satisfying G(x) = x/(x) > 0 could be used in order to 
maximize the dissipation. 


Then, we have 


which yields 


[F n \ = Q[(p]H = 2e r H and -Q[V r 0] = -V r e r 


[(ffj, + °e)\ ■ n ~ [Qp u \ = F mb = (Yo ~ e r )Hn - V r e r . 


( 20 ) 

( 21 ) 


and 


[£ 

Pi 


~\ n '( a V + a e) n + [ e es +^k B TCi J + Q5 n 0 -^P.-J-Ci 


+ 


\\u n \ 2 

QpO-n 

2 

P 


^(n-(a v )-n + 0-^p c ^-c^j 


+ 


Qp + 2Q p v n 

_Qp 

2 p 2 

P 


f Qp 

<p 
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1 ( ( \ V dip \] \ Q P ( 22 ) 

= -\n (o„) n + 0-^p c — Ci jj- [^j. 

Here we used the fact that v n is continous across the membrane. 

Model for fluid with variable density . To summarize, we have the following model for fluid with 
variable density passing through a deformable membrane 


-V • (sosfVcp 1 ) = ^ z^Cf, 


P ~ P( C 1 > C 2 ’ C w) 


¥ + F ' ( u±c « ± ) = "' [ D ‘ ( Vc ‘ + 
, dvr- , , . 

p± —— = v • (o-±) + v • o±), 


+ p ± V • ^ — 0, 


, in ^ 


in fi ± 


in D. ± 


in fi ± 


where viscosity shear stress is 


<r± = 2r7±D± + AHV ■ m±)/ - ^ k B Tcf + p± ^1 - ^ ^/ 
and Maxwell stress = e o ef : (V0 :!: ® — 1| V0| 2 /). with interface conditions on T 


-SoS^Vc/) 1 ■ n — C m [(/)], 


jf n + Cf (u± - v) ■ n = ^y{[k B T (In(q/c 0 ) + 1) + z^0 - p c ), 
Qp = -p 2 ff(*)j ^n- • n+ 0 -'YjPc-yC^J - ^ Q 2 p j 


Q P = ~P 


p ± (u ± -v)-n = Qr 


K + a e ] ■ n—[Q p u] = (y 0 - e r )Hn - V r e r 


U~ T — V T 


and some appropriate boundary conditions on dfL 
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Remark. Here it is easy to check that the model (23) with interface conditions (24) is Galilean 
invariant. For the transmembrane mass flux, note that (23) is a quadratic equation of Q p , which 
normally involves two solutions. If we denote 6 = [l/2p 2 ] « 1, one solution involves Qp as a small 
perturbation 


Q p = —p l K(x) 


+5p 2 K (x) j p 2 K (x) 


1 

P 



+ 0(5 2 ), 


which means Q p is mealy driven by the hydro and osmotic pressure difference. 
The other candidate solution makes Q p a large quantity 


Q P 


i i 

S p 2 K(x) 


p 2 K(x ) 



+ 0(5). 


which means Q p is meanly elriven by the density difference caul does not depend on the hydrostatic 
and osmotic pressure. The large Q p solution is not reasonable for a semi-permeable membrane flux 
because the mass flux increases as density difference decreases (8 -> 0). So the first solution will be 
chosen for our semi-permeable membrane mass flux. 


Remark Here we considered a simple membrane mechanical energy, only including the surface 
tension. The derivation could be generalized to more complicated membrane energy case, for 
example the Helfrich bending energy* 7 114 caul the neo-Hookean hyperelastic energy 115 , allowing 
models of fascinating biological phenomena like vesicle fusion 116 ' 118 to include the effects of ionic 
composition, membrane potential, and hydrostatic pressure that are known to have significant effects 
in experiment caul life. Details are in the Appendix B. 
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Permeable and Inflexible Membranes 


In this section, we assume that the membrane is inflexible v — 0 and the total energy and 
dissipation functional are same as in Eqs. (3) and (4). The flux interface conditions are reduced 
to p ± u ± = Q p n and )f ■ n + CfU^- ■ n — /j(x), on T . By using the similar calculation with the 
reduced interface conditions , Eq. (16) is reduced to 


dE tot 

dt 



<r^: dx 


^ J (a e -£(,£±(170* ® 70± -^|E0±| 2 /)j: Vudx 


*S LI (ypf) ■Jt dx - ll [fiil/idS 


(25) 


+ - a± i 


+ 




+ 




QpdS + I b . 


V■ u ± dx 


Comparing with the dissipation functional (4), we obtain the same results as in Eq. (19) without the 
last two equations for force balance on the membrane. 


Model for an inflexible semi-permeable membrane . To summarize, when the density of solution is 
a function of ion concentration, we obtained the following model for fluid passing through an inflexible 
semi-permeable membrane 


-17 • (£ O £±[70±) = ^ Zi ecf, 


dcr 


'-+V- (it±cf) = 7-1 £>f ( Vcf + t70± --^—cfV (p± 

dt K 1 J \ V k B T k B T V P , 

p±d ^r =v ' +17 ■ ( e o^r0± ® 70± - | E0±i 2 ), 


dp 1 * 1 

dt 


-1- • u* = 0, 


P — P( C 1 ’ C 2 ' Cv) 


in 


in 12 * 

(26) 


in 12* 

in fl ± 
in fl ± 
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with interface conditions on T 


-SoE^Vcf) 1 n = c m [(p\, 


jf n + Cfu 1 ■ n = ^ ( k B T (ln(q/c 0 ) + 1) + - p c ) , 


Q p - -p 2 K(x) 


-10 + 


srdip 

n '^' n “Zv 


PcC 


r 11 


p~u— = Q p n, 


and some appropriate boundary conditions on dfL 


(27) 


Remark Comparing with the results for flexible membrane in Eqs. (23)-(24), we find the only 
difference is that there is no force balance equation for the inflexible membrane. The inflexible 
membrane has an unknown external force that maintains the position of membrane. Only the Dirichlet 
interface condition is used to describe the solvent velocity passing through the inflexible membrane. 
In the case that the membrane is non-permeable for fluid, an equation of osmotic velocity is not needed. 
In this case, the membrane moves with the fluid, i.e., Q () — 0 (i.e. = v); only the force balance 

equation \a p + o e \ ■ n — (y 0 — e r )Hn — V r e r on the membrane is needed when the membrane is 
non-permeable for fluid. 


Model with Constant Solution Density 

For the normal biological problem, where the bulk ion concentration is in the range of 
1~300 mM, the density of solution has small variation. Concentrations are very much larger in and 
near ion channels, charged lipid membranes, binding proteins, enzyme active sites, and nucleic acids. 
In those situations solutions are also quite nonideal and so a different treatment is needed. 

In the bulk, we can take the zero order approximation, i.e. p ± — p°, where p° is a constant. 
Then the fourth equation in (1) is reduced to the imcompressiblility condition V • u = 0. The kinetic 

energy functional is reduced to be E kin — £+ J Q± ^p°\u ± \ 2 dx. The second term in dissipation 

functional vanishes. The interface conditions are reduced to 


r 



on r < 


D± n = Q{x), 
u ± • n — v ■ n — U, 
if n + Cf (m± -v) n = Ji(x), 
\^r) + ®"e] ' ^ — Fmb 


M ± •T = V •T 
V 


(28) 
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Here the main difference between Eq. (2) and Eq. (29) is that the velocity on the membrane is 

Qo 

continuous since the density is continuous. Here we replaced by new variable U. 


Permeable and Flexible Membranes 

We also first assume that the membrane is flexible. Then by using a similar calculation with 
interface conditions (28) and incompressibility, Eq. (16) is reduced to 


dE 


tot 


dt 


Vu ± dx 


H70±| 


-IL* 

i 

- X i n+ (** ~ e ° er± ( 70± 0v<pz 
+ X /+X VMt± ■ ^ dx - j X ds 

± n± i r i 

+ j (Fmb + [Q v 0 ]) • V TdS 

+ J ( Fmb • n + 1A] _ ( e r + Yo)H)V n dS 
+ j [F p ]udS + j [<p]^(C m [<p]-Q)dS + I b 


-I : Vu^dx 


(29) 


4 — hb + 4 b> chemical potential p. t reduces to /q; F n — Q4>H is the same as in Eq. (17); and 
F p = n (a v + a e ) n + ^e es + ^ k B T c f + Q d n (/}j. 


(30) 


Model for a deformable membrane . To summarize, the model for fluid passing through a deformable 
membrane is 


-V ■ (soefVc/) 1 ) = ^ Ziecf, 


dcr 


dt 


l -+V- (u±c±) = V ■ (D?(Vc? + ^c±l70±)), 


k R T 


in fl - 


in H ± 


du- 


£o£r 


p-^+Vp ± = E-(2?7 ± D±) + l7-(£ o £±E0 ± (g)l70 ± -^^|l70 ± | / /), in Q- 
dt 1 2 

V • u- = 0, 


in S2~- 


(31) 
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where the viscous stress for incompressible fluid is 

<r± = 217 ±D± - p±/ f 

and pressure p is introduced as a Lagrange multiplier of imcompressiblity. 
By taking p = p°, the interface boundary conditions on F are 

f dX 



-£ 0 £r ' 7 0 ± n= c m [0], 

ji n + cfiu 1 - v) ■ n = ([ fe B r (ln(Ci/c 0 ) + 1) + z £ e0]) f 

u ± —17 = Un, 

U — —K (x) [0 + n • o-jj • n], 

K + <Te\ ■ n = (.Yo ~ e r )Hn - V r e r , 

V 


(32) 


and some appropriate boundary conditions on dfl describing the experimental setup in a reasonably 
realistic, but idealized way. Note that the system (31) with interface condition (32) are same as the 
previous results 106 . 


Remark Other formulations could be used for the membrane of course. Even the tradtional Goldman— 
Hodgkin-Katz flux 106 

Ji — Pif(zie[(t)]/k B T)c[(e^ /kBT — 1) with /(x) = —^ > 0 ,for any x, 


could be used to replace the Hodgkin-Huxley flux Ji — 9l [p;]. 

v z i e ) 

What is recdly needed of course is a physical model of the water, mass, and electrical flux through the 
ensemble of channels, starting with atomic detail. That is not yet available, particularly because the 
macroscopic representation (of the ensemble of channels) needed here (i.e., in eq. (32)) should include 
time dependent gating phenomena that are not properties of a single permanently open channel. The 
Goldman-Hodgkin-Katz and Poisson Nernst Planck formulations are representations of the behavior 
of single permanently open channels. Those representations need to be supplemented by an explicit 
theory of time (and agonist or voltage) dependent opening and closing of both single channels and 
ensembles of single channels before they can sensibly be used in a model like eq. (32). Note two types 
of gating are involved, the stochastic open and closed gating of single channels, and the deterministic 
time dependent conductance ccdled gating in the Hodgkin Huxley formulation of ensemble properties. 
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Permeable and Inflexible Membranes 

In this section, we further assume that the membrane is inflexible, i.e. v = 0 . Similarly, the different 
between flexible and inflexible membrane is that we do not need the force balance equation on the 
interface. Then we could obtain the following model for fluid passing through a inflexible membrane 

in fl ± 

inO 1 (33) 

£r)£f- | I 

^|170±| 2 /), inH ± 

^ V ■ u— = 0, infl ± 

where the viscous stress for incompressible fluid is 

<r± = 2p±D± - p±/ 

with interface boundary conditions on T 


-V • (£ O £±[70±) = ^ Zi ecf, 


dcr 


z,e 


d l r +\ 7- (u±c t ±) = 17- (D?(Vc} + ^±170±)), 


du ± 


p “ + Vp ± — V ■ ( 2pD + V ■ ® 170* — 


-£ o £r' 7 0 ± - n = C m [0], 

jf n + CfM 1 • n = -^—([k B T(\n(c i /c 0 ) + 1) + Zje0]), 
m- - Un, 

U — — K(x)[0 + n ■ a v ■ n], 


(34) 


and some appropriate boundary conditions on dfl. 
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Historical Comments 


Osmosis occurs in so many contexts that no one, certainly not us, can grasp its totality. The 
temptation then is to discuss only what we have grasped and say nothing about what is beyond our 
reach. Because we do know something about some important applications—and think something is 
better than nothing—we write this section. We hope workers in different fields will learn of each 
other’s work, so they can benefit from each other’s knowledge, experience, and enthusiasm about 
osmosis. 

It is important also to understand that at least in biology many of the most important issues in 
osmosis remain open. We believe, without proof, that some of these issues will be easier to resolve 
when ionic solutions are treated as complex fluids with mathematics that is consistent, with all 
variables satisfying all field equations and interface boundary conditions with one set of unchanging 
parameters. We understand, however, that no one can know the consequences of a consistent analysis 
until it is actually performed in the experimental, technological and biological systems of interest. 
Details matter! 

An enormous amount of work is being done on the theory of technological osmosis, i.e., 
mathematics of fluidics. Our knowledge of fluidics is too modest to allow thoughtful citation here: 
fortunately these applications are easy to locate with literature search methods. A search today (June 
2, 2018) found two significant new papers just this week. 

The classical biological literature is older, and in danger of being lost to modem generations, 
particularly those who know more of Navier Stokes and fluidics than of the kidney or epithelia. So we 
provide key references to interesting biological applications. It is important to understand that these 
areas are central to a wide range of biological research. They are not isolated special cases. 

Much of classical physiology, described in detail in textbooks 21 " 24 concerns organs that depend on 
osmosis. Mori 106 ' 119 has started a consistent analysis of osmosis and our work should be viewed as an 
extension of his. 

The reviews of Boulpaep 120 , AE Hill 121 , and Pohl 122 provide good entries to the field. Hill is 
particularly useful for showing the substantial controversies and their history. The analysis of osmosis 
in the lens, mostly from the laboratory of Richard Mathias, is notable for its success, combining field 
theory, molecular biology, structural biology, and measurements of hydrostatic pressure, electrical 
impedance, and fluxes to provide a coherent view of how the lens uses osmosis to stay alive. 123 " 125 
That approach, expanded perhaps to explicitly use the theory of complex fluids, will help to resolve 
the many controversies we believe. 

A great deal of work has been done on fusion of vesicles to membranes, because of its wide 
biological role (e.g., in the ~10 14 chemical synapses in our brains, with thousands of vesicles in each 
synapse) and medical importance in the entry of viruses into cells. 116,117, l26 ‘ 127 Most of that work only 
considers the elastic properties of membranes. We suggest that energy sources like diffusion, electrical 
potential, and convection should be included in the analysis of vesicle fusion—and membrane flow in 
general—to see if they are used by biology for its purposes. Experiments suggest they are. 

Modem work in molecular biology and water flow is focused on the aquaporins that are thought 
to provide channels for water flow 128,129 , or sensors controlling water flow 130,131 . Some work seeks to 
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study osmosis in individual channels, but comparison with similar work studying current flow through 
channels suggests that much remains to be done. In our view, the full range of forces—electrical, 
diffusional, convective, steric, and chemical—are likely to be involved in flow through water channels 
and so must be included in models and simulations, at least to begin with, until we discover the 
variables that evolution uses to control water flow on the atomic scale in the channels of aquaporins. 

Sadly, osmosis and fluid flow in plants is not an area we know enough about to cite intelligently. 
Osmosis/fluidics reaches its most impressive heights in trees (nearly 100 meters in redwoods). Water 
flow in the xylem of trees—from roots to leaves in the crown—certainly needs analysis with the theory 
of complex fluids, if that has not already been done. There are also innumerable examples of 
osmosis/fluidics we do not know enough about to even name. All the more reason for readers interested 
in complex fluids to learn their names and study how they work. 

In each case, it seems clear to us that treating osmosis as the movement of a complex fluid is 
likely to be useful, even necessary, in biology, and everywhere else. 
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Appendix A: Proof of Lemma 


Proof. By using the Batchelor’s method 115,132 , we have 

— ( n)dS + — ( dS)n = — ( ndS) — (V ■ v)ndS — ( Vv ) • ndS 
dt dt dt 

= (l V vn— ( d n v ■ ri)n)dS — V r i7 • ndS = (V r • v)ndS — V r v • nds 

and —dS = V r -vdS. 

dt 1 

Then we have — - — V r v ■ n. 

dt 1 

d n 

Finally, by the definition of —, we have 
d n n dn 

~dF = ~dt~^ Vr ' Vr ^ n = ~ VfV ' n ~( Vr ' V ^ n 

= —V r (y n n + v T )-n - V r n • v r = -V r v n - V r (n • v T ) = -V r v n . 

Here we used the symmetry of curvature tensor V r n and \n\ 2 = 1. 
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Appendix B: Membrane mechanical properties 

Membrane mechanical properties can be included in the analysis by modify Eq. (7). Assuming 
that E m is the energy induced by mechanical properties of the membrane, the fourth term of Eq. (7) 
/ 4 yields 


dE m 

dt 


4 - T- f (.Yo + C m [(l)] 2 )dS + 
at j r(t) 

= | C m [0]^([0])d5-J ((e r Hn — F m )) • vdS. 


(SI) 


where we used the results 87 = f . • vdS and denotes F m — ^ rrL 

dt J r(t ) sr m sr 

Then the rate of total energy change is 


dE 


tot 


dt 




-Z L* 




+ X /+Z v ^’ ; ^ dx ~ / Z^^ ds 


+ 12 


|E0±| 


■/ : Vu^dx 


+ 




dx 


+ | (Fmb -F m + [QV0]) • v T dS 
+ f (Fmb ■ n + [F n ] - (e r H - F m ■ n))v n dS 

+ {[^] <3pdS + / [^ (Crn[0]_Q)d5_4 
Then we obtain the force balance equation on a membrane with mechanical properties. 


K + a e ] ■ n - [Q p u\ = -e r Hn - V r e r + F m . 


(52) 


(S3) 
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